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Abstract 

Several devices exhibiting memory effects have shown up in nonlinear circuit theory 
in recent years. Among others, these circuit elements include Chua's memristors, as 
well as memcapacitors and meminductors. These and other related devices seem to 
be beyond the, say, classical scope of circuit theory, which is formulated in terms of 
resistors, capacitors, inductors, and voltage and current sources. We explore in this 
paper the potential extent of nonlinear circuit theory by classifying such mem-devices 
in terms of the variables involved in their constitutive relations and the notions of the 
differential- and the state-order of a device. Within this framework, the frontier of 
first order circuit theory is defined by so-called hybrid memristors, which are proposed 
here to accommodate a characteristic relating all four fundamental circuit variables. 
Devices with differential order two and mem-systems are discussed in less detail. We 
allow for fully nonlinear characteristics in all circuit elements, arriving at a rather ex- 
haustive taxonomy of C 1 -devices. Additionally, we extend the notion of a topologically 
degenerate configuration to circuits with memcapacitors, meminductors and all types 
of memristors, and characterize the differential-algebraic index of nodal models of such 
circuits. 
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1 Introduction 



Broadly speaking, nonlinear circuit theory is concerned with the study of constrained ordi- 
nary differential equations involving time and four m-dimensional variables q, tp, i, v (stand- 
ing for charge, flux, current and voltage, respectively), with the following restrictions: 

(a) the 2m differential relations q' = i, ip' = v always hold; 

(b) the vectors v, i, satisfy a total amount of m linearly independent relations Bv = 0, 
Di = coming from Kirchhoff laws; 

(c) the characteristics of devices define m additional relations among the circuit variables. 
This means that nonlinear circuit models can be generally written in the form 

q = % (la) 
<p' = v (lb) 
= f(q,<P,i,v,t), (lc) 

where / captures Kirchhoff laws but also the constitutive relations of all circuit devices. 
Aside from the different circuit topologies, reflected in the form of the loop and cutset 
matrices B, D arising in item (b), the differences between circuit families come from the 
devices' characteristics referred to in (c). 

In classical circuit theory, the devices' characteristics just involve two out of the four 
variables mentioned above in one of three ways, relating current and voltage in resistors 
and controlled sources, charge and voltage in capacitors, and flux and current in inductors. 
In general, these relations may be nonlinear and/or involve time explicitly; specifically, the 
current or the voltage in independent sources is an explicit function of time. This provides a 
setting which models a great variety of circuits arising in electrical and electronic engineering; 
in particular, many semiconductor devices, including transistors, may be accommodated in 
this framework by means of equivalent circuits including controlled sources. 

Memristors. The memory-resistor or memristor is changing this picture substantially. 
This device, whose existence was predicted by Leon Chua in 1971 [S] and which actually 
appeared at the nanometer scale in 2008 [IS], is defined by a nonlinear charge-flux rela- 
tion. The potential applications of this device in the design of non-volatile memories, signal 
processing, adaptive and learning systems, reconfigurable nanoelectronics, etc., might make 
the memristor and related devices play a very significant role in electronics in the near 
future, specially at the nanometer scale. A lot of research is focused on this topic; cf. 

[210[I21[lIl[23[283EniE 

Not only regarding applications, but also from a theoretical point of view this scenario 
poses challenging problems. Several extensions are being developed; some remarkable ones 
are related to the memristive systems of Chua and Kang [TO] and the mem-devices (mem- 
capacitors and meminductors) recently introduced by Di Ventra et al. [IT] . Somehow, the 
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fundamentals of nonlinear circuit theory are affected, since its concepts and models seem to 
be moving beyond their classical limits, and several questions arise. One may wonder about 
the limits to which this framework may extend, but also how to accommodate these new 
devices in a comprehensive taxonomy, and which aspects should articulate such a taxonomy. 
In the present paper we will try to answer some of these questions. 

In this regard we will find it useful to describe the memristor (say, in a charge-controlled 
setting) not by means of its flux-charge relation ip = <p(q) but, instead, via the differentiated 
relation <p' = 4>'(q)q' , that is, 

v = M{q)i, (2) 

where M(q) = (j)'{q) is the so-called memristance. The reason to do so is that, in the 
description of the dynamics of memristive circuits, there is no reason to keep track of both 
the flux and the charge as dynamic variables since their values are constrained by the relation 
ip = <p(q). The memristor has order (or, to avoid terminological ambiguities, state order) 
one. This means that every memristor introduces one, but not two, degree(s) of dynamic 
freedom, associated with the charge q; describing the device characteristic via (T5]) we get rid 
of the flux ip. A detailed discussion can be found in subsection 12.11 

Now, the key remark in our discussion is to look at (j2]) just as a relation involving three 
out of the four variables q, <p, i, v (all but <p), regardless of its specific form; actually, as 
in [JT] we will work with so-called q-memristors, defined (in a time-invariant setting) by a 
general, fully nonlinear characteristic of the form v = rj(q,i). 

The corresponding relation for y9-memristors is % = (((f,v) and, in particular, the char- 
acteristic of Chua's flux-controlled memristors [8] reads as 

i = W{cp)v, (3) 

in terms of the memductance W(ip); the identity ([3]) again involves three out of the four 
variables mentioned above (in this case, all but q). 

Memcapacitors and meminductors. In a natural way, this makes us wonder about the 
remaining relations involving the other two combinations of three variables, namely, q, ip, v 
and q, <p, i. This will lead to the (voltage-controlled) memcapacitor, with a characteristic of 
the form 

q = C m (<p)v, (4) 

and the (current-controlled) meminductor, for which, 

if = L m (q)i. (5) 

These devices have been introduced by Di Ventra et al. in [11] and are discussed in subsection 
I2.2[ using again more general, fully nonlinear characteristics. 
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Hybrid memristors. Finally, for the sake of completeness one may think about a char- 
acteristic actually involving all four variables q, <p, i, v. For reasons detailed later, we will 
consider two different settings, defined by the fully nonlinear characteristics v = ip(q,ip,i) 
and i = £(q,(p,v). In particular, when these relations are linear in i and v, respectively, we 
get 

v = M h (q,p)i (6) 

and 

i = W h (q, <p)v. (7) 

We will call these circuit elements hybrid memristors, since their memristance and memduc- 
tance involve both the charge q and the flux (p. These novel circuit elements, together with 
their fully nonlinear variants, are introduced here for mathematical completeness, but they 
also provide a natural framework to accommodate actual physical devices in which different 
memory effects (namely, memristive, memcapacitive and/or meminductive features) coexist. 
Find details in subsections 12.31 and 12.41 

These relations comprise all possible combinations of the four variables q, <p, i, v. This 
means that the aforementioned devices and their time-varying analogs exhaust what may be 
called first order circuit theory, within the setting defined by items (a), (b) and (c) on p. [2j 
Here the term "order" is used to mean differential order (cf. Section [2]) and, accordingly, the 
expression "first order circuit" reflects the fact that only the branch currents and voltages i, 
v and their first integrals q, p are involved (compare Definitions [1] and [8] in Sections [2] and 
[3]). Certainly, the reader should not misunderstand our use of this expression with that in 
elementary circuit theory, referring to a circuit with only one reactive element. 

Beyond this framework, higher order devices involve additional variables such as the first 
integrals of the charge and the flux (second integrals of the branch current and voltage), to be 
denoted by a and p, respectively. With the use of these additional variables we may accom- 
modate devices such as charge-controlled memcapacitors and flux-controlled meminductors 
Another extension of the theory is related to the mem-systems introduced by Chua 
and Kang [10] . Roughly speaking, the restriction on the form of the differential relations 
considered in (a) is no longer assumed in these systems; by contrast, new variables without a 
physical meaning in classical circuit-theoretic terms are allowed, providing dynamical states 
on which the resistance, capacitance and inductance (or better, the memristance, memca- 
pacitance and meminductance) may depend. Higher order devices and mem-systems are 
briefly discussed in Section [3j a deeper analysis of them is in the scope of future research. 

The use of these types of devices in applications, as well as the numerical simulation of 
the dynamics of nonlinear circuits including them, requires some effort also at the circuit 
modelling level. We will consider circuit modelling aspects in Section |H Specifically, we 
analyze there the differential-algebraic models arising from the nodal analysis of nonlinear 
circuits including all possible types of first order devices. Our goal in this regard will be the 
characterization of the index of these models, a notion which defines many key properties 
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of differential-algebraic systems jU [HJ CEO [231 EH ES]- This way we will arrive at a general 
index characterization (cf. Theorems [1] and [2] in Section H|) from which the results discussed 
in [121 EH SSI S3 HE] can be derived as particular cases which apply to circuits including only 
certain types of devices with differential order one; another particular case of importance 
concerns circuits with memcapacitors and meminductors, whose index is so far unexplored 
in the literature. We will also extend the notion of a topologically degenerate configuration 
to this broader setting. 

2 First order devices 

Definition 1. A circuit device is said to have differential order one if it is defined by a 
C 1 - characteristic of the form 

h(q,<p,i,v,t)=0, (8) 
where at least one of the partial derivatives h q , h v does not vanish identically. 

Note that both h q and may of course vanish for specific values of the circuit variables. 
Definition [1] is supported on the fact that when both h q and vanish identically, the only 
branch variables involved are the current and/or the voltage. This amounts to what may 
be called (differential) order zero devices, not involving any dynamics. These are nonlinear 
resistors, and voltage and current sources. By a nonlinear resistor we mean any device 
relating current and voltage in an algebraic (i.e., non-differential) manner, such as a diode. 

The use of the term "order" is worth a digression. This term has different senses in 
mathematics, and at least two apply in our context. On the one hand it means the order of 
derivation of a system of differential equations; in this sense, Newton's law yields a second 
order system. We will use the expression differential order to mean this. On the other hand, 
"order" is also used to mean the number of dynamic variables in a system of differential 
equations; the term is often used in this sense in circuit theory, to refer to the number of 
state variables associated with a given device. We will describe the latter as the state order 
of a device. Throughout the document and when no label is used, by "order" we refer to the 
differential order; e.g. a first order device is a device with differential order one. 

As detailed in the sequel, the list of devices with differential order one will include 

• capacitors, inductors, g-memristors and 92-memristors, all of them with state order one; 

• voltage-controlled memcapacitors, current-controlled meminductors, and hybrid mem- 
ristors, with state order two. 

In practice, at least two of the derivatives h q , h^, hi, h v do not vanish identically in first order 
devices. It is also worth noting that in Definition [T] we implicitly assume the arguments of h 
to be one- dimensional; however, allowing q, ip, i and v to take vector values, this definition 
easily accommodates coupled devices. In this case the non-vanishing requirement on h q 
and/or h v must be replaced by the non-singularity of the corresponding matrices of partial 
derivatives. This remark will apply throughout the paper, often without explicit mention. 
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2.1 q- and (^-memristors 

The electromagnetic relations q' = i and ip' = v link the pairs of fundamental circuit variables 
q — % and p — v. First order devices involving the other combinations of two out of the four 
variables q, p, i, v are the capacitor, relating q and v, the inductor, which involves <p and i, 
and Chua's memristor, whose characteristic relates q and p. 

The complete dynamical description of the first two devices requires including the dif- 
ferential relation q' = i for the capacitor (note that ip does not appear in its constitutive 
relation), and ip' — v for the inductor (q not being involved). Regarding the memristor, 
even though both q and p arise in its characteristic, these variables are linked together and 
therefore only one of them introduces a dynamical degree of freedom in the circuit. For this 
reason, it will be preferred to get rid of either the flux or the charge by means of a relation 
formulated in terms of the other three variables, as detailed below. This way, not only the 
capacitor and the inductor but also q- and <^-memristors will have state order one. 

Definition 2. A g-memristor is a device with differential order one, governed by the relations 

q = i (9a) 
v = r){q,i,t), (9b) 

where rj is a C x -map for which neither of the derivatives 7] q , rji vanishes identically. 

When 7] is time-invariant and linear in i, we get Chua's characteristic v = M(q)i, where 
M(q) is the memristance [8]. The identity q(t) = J_ oo i(r)dT shows that this volt age- current 
relation keeps track of the device history; for this reason Chua proposed the name memory- 
resistor, or memristor for short. In the literature, this device is said to be current- controlled 
but also charge-controlled, because of the form of the map p = <fr(q) whose time derivative 
yields the voltage-current relation v = M(q)i. 

For fully nonlinear memristors of the form (Q the incremental memristance is defined 
as the derivative r]i(q,i,t) [H] ; the device is called strictly locally passive if rji(q,i,t) > 
for all (q,i,t) (or if the matrix rji is positive definite in coupled cases). The requirement 
that rji does not vanish identically distinguishes the device from a nonlinear capacitor. In 
turn, the non- vanishing condition on rj q makes this device actually different from a nonlinear, 
current-controlled resistor. Note also that the fully nonlinear form (Q makes it possible to 
accommodate devices displaying memristive effects but whose characteristic does not arise 
as the time derivative of a p-q relation, contrary to Chua's memristor. 

Definition 3. A y?-memristor is a device with differential order one, governed by 

(ff = v (10a) 
i = ((p,v,t), (10b) 

where ( is a C l -map such that neither of the derivatives C>; Q v vanishes identically. 
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A time-invariant <y2-memristor for which (llObj) is linear in v (hence reading as % = W(p)v ) 
amounts to Chua's flux-controlled memristor [8j, W(<p) being the memductance. In general, 
the incremental memductance is the derivative ( v (<p, v, t), and the device is said to be strictly 
locally passive if this derivative is always positive. Again, the non-vanishing requirements 
on and ( v make this device different from a voltage-controlled resistor and an inductor, 
respectively 

Nonlinear, C l circuits composed of (current- and voltage-controlled) resistors, voltage and 
current sources, capacitors, inductors, and q- and <£>-memristors display the property that 
neither the differential nor the state order of their devices exceed one. Several dynamical 
properties of this type of circuits are discussed in [UJ. In the next subsections we discuss 
the main features of devices with differential order one and state order two. 

2.2 Memcapacitors and meminductors 

The remaining characteristics involving three out of the four variables q, ip, i, v naturally 
lead to the voltage-controlled memcapacitors and current-controlled meminductor discussed 
below. Both devices will have state order two. We consider fully nonlinear characteristics 
for memcapacitors and meminductors; when the maps u and 9 below are time-independent 
and linear in v and i, respectively, we will get the circuit elements introduced by Di Ventra 
et al. in [TT] . 

Definition 4. A voltage-controlled memcapacitor is a device with differential order one, 
defined by 

q = i (11a) 
<p' = v (lib) 
q = u(<p,v,t), (11c) 

where u is a C 1 -map and neither of the derivatives u v , u v vanishes identically. 

When u is time-invariant and linear in v, (II left reads as q = C m (ip)v, where C m is the 
memcapacitance [TT] . The distinct feature of this device is that the memcapacitance depends 
on the state variable (p(t) = J_ oo v(r)dT, so that the relation q{t) = C m ( f v(r)dr)v(t) 
reflects the device history. Be aware of the circuit-theoretic meaning of this state variable; 
when the memcapacitance is allowed to depend on abstract state variables we are led to 
the more general setting of memcapacitive systems (see [TT] and Section [3] below) . It is also 
worth mentioning that the relation q = C m ((p)v arises as the derivative of a characteristic 
a = fJ,(ip), where a is the time integral of q. This yields a 1 = q = (j,'((p)ip' = C m (<p)v. Notably, 
the device can be described without recourse to the (second order) variable a (cf. (|lip ). in 
contrast to the charge-controlled memcapacitors discussed in Section El 

For an arbitrary C 1 -map u, the incremental memcapacitance C m is defined as the deriva- 
tive u v and, in general, depends on (ip,v,t). Such a device need not come from a a — ip 
relation, and the requirement that neither u 9 nor uj v vanishes identically makes it actually 
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different from a capacitor or a memristor, respectively. An instance of a fully nonlinear 
memcapacitor arising in a Josephson junction model can be found in subsection 12.41 

Definition 5. A current-controlled meminductor is a device with differential order one, 
governed by 

i (12a) 

v (12b) 

9(q,i,t), (12c) 

where 8 is a C l -map for which neither of the derivatives 9 q , $i vanishes identically. 

Again, when the map 6 is linear in i and does not depend on t, we get the characteristic 
<p = L m (q)i = L m (J i{r)dr)i considered by Di Ventra et al. in [TT]. Such a characteristic can 
be obtained as the time derivative of a C 1 -relation p = n(q), where p is the time integral 
of ip; note however that the description (fT2"j) does not involve p (compare with the flux- 
controlled meminductors in Section [3]). Now L m is the meminductance. In general, the 
derivative 9i(q,i,t) is the incremental meminductance. The non-vanishing of this derivative 
makes the meminductor different from a memristor and, similarly, the assumption that the 
partial derivative 9 q does not vanish identically makes the device different from an inductor. 

The appearance of both q and (p in the characteristics flllcl) and f!12cj) imply that both 
variables must be present in the dynamical description of memcapacitors and meminduc- 
tors and, therefore, that these devices have state order two; this is an important difference 
with capacitors, inductors and q- and 92-memristors, for which one dynamic variable suffices 
(together with i and v) to describe the device behavior. 

2.3 Hybrid memristors 

In light of the characteristics ( 19bj) . fllObl) . flllcl) and (I12cj) . from a mathematical point of 
view it is somehow natural to complete the picture by considering a relation which involves 
all four variables q, ip, i and v. We present below two different settings (dual to each other) 
in which these four variables might actually arise. 

These (so-called hybrid) memristors may account for physical devices in which memory 
effects of different nature coexist. The simultaneous appearance of memristive, memcapaci- 
tive and/or meminductive phenomena has been discussed by different authors. For instance, 
the coexistence of memristive and memcapacitive effects has been reported to follow from the 
formation of local dipoles in nanoscale resistors [TT] , and may also arise in metal-insulator- 
metal thin films having thickness between the nanometer and the micrometer scales [29J. 
Capacitive and memristive effects coexist with the nonlinear inductive nature of a Joseph- 
son junction in accurate models of this device [2T]. From a modelling point of view, hybrid 
memristors also allow for simplified descriptions of combinations of more basic devices. Find 
examples in subsection 12.41 below. 



Q = 
V' = 

<P = 
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Definition 6. A current-controlled hybrid memristor is a device with differential order one, 
defined by the relations 

q' — i (13a) 
<p' = v (13b) 
v = ip(q,<f,i,t), (13c) 

where ip is a C 1 -map and none of the derivatives ip q , ip^, xpi vanishes identically. 

The terminology is motivated by the fact that, in cases in which ( I13c[) is time- independent 
and linear in i, this relation takes the form v = M^(q,if)i, providing an analog of Chua's 
memristor in which the so-called hybrid memristance Mh(q,(p) now depends on both the 
charge and the flux; both variables introduce memory effects on the device, because of 
the relations q(t) = J ^i^dr, <p(t) = /_ 00 i'(r)£Zr. In general, the incremental hybrid 
memristance is defined as the derivative ipi, and may depend not only on q and ip but also 
on the current i and on t. The hybrid memristor is said to be strictly locally passive when 
the incremental hybrid memristance is positive (or positive definite in coupled cases). 

Definition 7. A voltage-controlled hybrid memristor is a device with differential order one, 
governed by 

q = i (14a) 
cp' = v (14b) 
i = €(q,<P,v,t), (14c) 

where £ is a C^-map for which none of the derivatives £ g; £ v vanishes identically. 

Cases in which £ is linear in v and time-invariant yield i = Wh(q, (p)v, where is now 
the hybrid memductance. In this situation, when is non-singular for all values of q, ip, 
then obviously the device can be recast in a current-controlled form and vice-versa. Again, 
in general the incremental hybrid memductance is defined by the derivative £ v (q,ip,v,t) 
and, as before, the device is said to be strictly locally passive when the incremental hybrid 
memductance is positive (or positive definite). 

The non- vanishing requirements on the derivatives within the characteristics f)13c[) . (I14cl) 
distinguish these devices from the q- and yj-memristors, memcapacitors and meminductors 
discussed above. Actually, the non- vanishing of these derivatives implies that, at least locally, 
both the charge and the flux can be written in terms of the remaining circuit variables. 
Focusing the attention on the current-controlled case, the fact that ip q ^ makes it possible 
to recast fl!3cj) . via the implicit function theorem, as 

q = a((p,i,v,t), (15) 

and the device exhibits a (generalized) memcapacitance 

a v ((f, i, v, t) = 4>- l (a(ip, i, v, t), <p, i, t). 
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Similarly, the non-zero nature of the partial derivative ip^ yields 

<P = P(q,i,v,t), (16) 
for some locally defined map /3, the (generalized) meminductance being 

Pi(q,i,v,t) = -i)~ l (q,(5{q,i,v,t),i,t)i)i(q,l3{q,i,v,t),i,t), 

as a consequence of the implicit function theorem. Similar remarks apply to voltage- 
controlled hybrid memristors. 

Noteworthy, removing the non-vanishing requirements on the derivatives, the relations 
(I13c| . (j!4cj) . ( IT5|) and (|T5|) account for the constitutive relations of all previous devices. From 
a dynamical point of view, a remarkable difference between hybrid memristors and both q- 
and (/9-memristors is that hybrid ones have state order two, since both the charge and the 
flux are necessarily involved in the description of the device. This implies, for instance, that 
the dynamics of a circuit with just one hybrid memristor and without reactive elements lies 
on the plane. The actual dynamical properties of circuits with hybrid memristors are in the 
scope of future research. 

2.4 Examples 

Example 1. Our first example attempts to illustrate how the use of fully nonlinear charac- 
teristics may allow for simplified descriptions of certain devices. Specifically, we will provide 
an accurate description of a Josephson junction by means of a fully nonlinear memcapacitor, 
which accounts for all parasitic effects, connected in parallel to a nonlinear inductor (cf. 
Figure [TJ below). 




Figure 1: (a) Josephson junction, (b) Equivalent circuit with a fully nonlinear memcapacitor. 

As detailed in [21], realistic models of a Josephson junction should take into not only the 
usual nonlinear inductive relation 

ij = h sm(k ipi) 

for certain physical constants I , k (see e.g. JD]), but also the presence of memristive, resistive 
and capacitive effects; an accurate equivalent circuit of the Josephson junction is defined by 
the parallel connection of these four elements (cf. [2T]), as depicted in Figure [T](a). 
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The device on the left of Figure Ufa) is a yj-memristor of Chua type, which as reported 
in [21] captures the presence of a small current component given by 

% w = h cos(ki(p w )v, 

again for certain constants I±, fci; here v is the port voltage. 

A linear resistor and a linear capacitor in parallel are also present in the description 
provided in [21J, being denoted by G, C, respectively, in Figure [fla). These elements are 
defined by the relations i g = Gv and q c = Cv. 

Now, the parallel connection of the memristor and the resistor is obviously governed 
by the current-voltage relation i wg = I\ cos(kiip w )v + Gv or, equivalently, by a charge-flux 
characteristic of the form 

Qwg = (h/h) Bm(kiip w ) + G<f w , (17) 

where we use the fact that the memristor flux is the time-integral of the port voltage v. 
The expression depicted in (|T7|) shows that the parallel connection of the memristor and the 
resistor is itself a <^-memristor. 

In turn, the parallel connection of the original memristor, the resistor and the capacitor 
can be described as a single device by setting q = q wg + q c , ip = ip w . Indeed, denoting by i 
the sum of the currents through the memristor, the resistor and the capacitor, we get 

q' = i (18a) 
up' = v (18b) 
q = (I x /kx) sm(kiLp) + Guj + Cv. (18c) 

This corresponds to a time-invariant, voltage-controlled memcapacitor for which the con- 
stitutive relation q = u(ip,v) in (II left takes the specific form depicted in (I18cj) . The corre- 
sponding equivalent circuit for the Josephson junction is displayed in Figure QJb). 

Example 2. Simple instances of hybrid memristors are defined by the connection of a 
charge-controlled and a flux-controlled memristor of Chua type; cf. Figure [2], where the 
charge-controlled memristor and the flux-controlled one are painted in green and yellow, 
respectively. In spite of their simplicity, these examples are not trivial: the goal is, again, 
to provide a dynamical description of each connection as a two-terminal device in terms of 
a single set of variables q, ip, i, v. 

In both cases, the subscripts 1 and 2 will correspond to variables associated with the 
charge- and the flux-controlled memristor, respectively. The charge-controlled memristor is 
assumed to be governed by a relation of the form ipi = 4>(qi), with memristance M(qi), and 
the flux-controlled one is defined by q<i = 7(^2), with memductance W(<p2)- Within Figure 
(resp. 12(b)), the memductance W{uo2) (resp. the memristance M(q\)) is assumed not 
to vanish. 

In the series connection of Figure EJ^a), elementary circuit theory yields % = %\ = i%, 
v = Vi + V2- In order to arrive at a dynamical description in terms of i, v and single variables 
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Figure 2: Examples of (a) a current-controlled and (b) a voltage-controlled hybrid memristor. 



q, (p, we set q = qi and (p — <px + <^ 2 . Obviously, this yields the relations q' = i and (p' — v 
but, more important, allows for the description of the voltage-current relation in the form 

v = M( gi )i + (W(<p 2 ))-H = [M(q) + (W(<p - y?!))- 1 ]* = [M(q) + (W(cp - ^(gi)))" 1 ]^ = 
= [M(q) + (W^- < P(q)))- 1 ]t. 

This corresponds to a current-controlled hybrid memristor for which the characteristic is 
linear in i, the hybrid memristance being 

M h {q^) = M{q) + {W^-cP(q)))-\ 

The physical meaning of the variable q is worth some additional remarks. Because of the 
identities q[ = i\ = i% = q' 2 , the charges q% and q 2 differ in a constant which would be fixed 
by the initial conditions in Chua's memristors. In turn, q is defined up to a constant, and 
therefore can be understood to describe any of both charges except for a fixed quantity. 
Mathematically, setting q = q\ + k for any real constant k, we would get a dynamical 
description of the device which, except for an affine change of coordinates, amounts to the 
one above; the description above assumes k = so that q actually equals the charge q\. 

It is also worth noting that if M(q) + (W(ip — does not vanish for all values 

of q, <p, the device also admits a voltage-controlled description. This would be the case, in 
particular, if the original Chua memristors are strictly locally passive (i.e. if M > 0, W > 
everywhere); in this case the hybrid memristor would itself be strictly locally passive since 
Mh would be strictly positive. 

Setting q = qi + q 2 , <p = y?2, the reader can proceed analogously in order to describe 
the parallel configuration of Figure 12(b) as a voltage-controlled hybrid memristor, with a 
characteristic linear in v and hybrid memductance 

W h {q, if) = {M(q - 7(^))) _1 + W[tp). 

These examples show that, even in simple cases, these devices pose interesting problems 
from the modelling point of view. Needless to say, the scope of hybrid memristors goes 
however beyond these academic examples. The definition of this circuit element is aimed at 
modelling devices which are not reducible to a simple connection of q- and ^-memristors, 
and which may capture the coexistence of different memory effects. In such devices the char- 
acteristics need not be linear in i and v , and the charge and the flux within the (incremental) 
memristance or memductance might interact in more intricate ways. 
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3 Higher order devices and mem- systems 

The devices discussed in Section [2] virtually fill the scope of circuit theory within the limits 
defined by the use of the variables q, ip, i, v and under the restrictions stated in items 
(a), (b) and (c) on p. [2J The frontier of nonlinear circuit theory in this setting is defined 
by hybrid memristors, which involve all four fundamental circuit variables. However, other 
devices recently introduced are located beyond these limits: on the one hand, certain circuit 
elements involve not only those four variables but also the time-integrals a, p of the charge 
and the flux, respectively; these include charge-controlled memcapacitors and flux-controlled 
meminductors [UJ, but also devices directly relating a and p, as proposed in [3J; all of them 
would have differential order two. On the other hand, allowing for new state variables, 
without the restriction imposed by (a), leads to the so-called mem-systems, which in general 
cannot be accommodated in the framework of Section [2J In this Section we briefly discuss 
some modelling aspect of these devices, leaving a deeper analysis for future research. 

3.1 Second order devices 

Definition 8. A circuit device is said to have differential order two if it is defined by a 
C l - characteristic of the form 

H <T >p>q> ( p>h v > t ) = Q > ( 19 ) 

where at least one of the partial derivatives h a , h p does not vanish identically. 

The non- vanishing of h a and/or h p implies that at least one of the differential relations 
a' = q (and in turn q' — i) or p' — ip (together with <p' = v) must be included in the 
dynamical description of the device. The differential order two nature of them stems from 
the identities a" = i and p" — v. Akin to differential order one devices, in practice at least 
two of the derivatives h a , h p , h q , h^, hi, h v will not vanish identically. 

Definition 9. A charge-controlled memcapacitor is a device with differential order two, 
defined by the relations 

a' = q (20a) 
q' = % (20b) 
v = u(a,q,t), (20c) 

where v is a C l -map for which neither of the derivatives v a , v q vanishes identically. 

Definition 10. A flux-controlled meminductor is a device with differential order two, gov- 
erned by 

p' = tp (21a) 
p = v (21b) 
i = X(p,f,t) (21c) 
where x is a C x -map such that neither of the derivatives x P , X<p vanishes identically. 
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From the descriptions (|20|) and (I2ip it follows that charge-controlled memcapacitors and 
flux-controlled meminductors have state order two; note that neither p nor (p (resp. er, q) 
arise in the characteristic f |20cj) (resp. (121 cj) ) . 

In particular, when z/ and \ are time-invariant and linear in q and respectively, one 
gets the devices introduced by Di Ventra et al. in for which the characteristics (I20cl) 
and (I21c|) read as 

v = E m (a)q (22) 

and 

% = K m {p)ip. (23) 

Here E m and TZ m are the inverse memcapacitance and the inverse meminductance, which 
introduce memory effects in the circuit because of their dependence on a and p, respectively. 

Noteworthy, the relations fl22l) and fl23l) arise as the differentiated form of certain map- 
pings ip = 7(0") and q = S(p), via the relations o' — q, p' — ip. In a natural way this 
leads to other second order devices, such as those relating a and p (cf. [3]). The analysis of 
the dynamics of nonlinear circuits including these or other higher order devices is an open 
problem, beyond the scope of the present paper. 



3.2 Mem-systems 

Mem-systems, originally introduced by Chua and Kang in [10J, are characterized by the 
removal of the classical electrical meaning of the state variable which introduces memory 
into the different devices. This means that new variables arise in the model, and that the 
specific form of the differential relations within item (a) on p. [2] is no longer assumed. For the 
sake of brevity, we restrict the discussion to the mem-systems which generalize the devices 
considered in Definitions [2j [3j H] and [5j using again fully nonlinear characteristics. 
Memristive systems are defined either by a system of the form 

x' = f(x,i) (24a) 

v = fj(x,i,t), (24b) 

or by 

y' = g(y,v) (25a) 

i = ((y,v,t). (25b) 

The distinct feature with respect to q- and yj-memristors is that now the memristance and 
the memductance depend on arbitrary state variables x, y, respectively, which do not even 
need to be one-dimensional; making x = q, f(q,i) — % and y = (p, g(<p, v) — v we get the q- 
and yj-memristors introduced in Definitions [2] and [31 Systems in which the characteristics 
f!24bl) and f!25bl) amount to v = M(x, %)% and i = W(y, v)v, respectively, describe the settings 
originally discussed by Chua and Kang in [10] . 
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Analogously, a memcapacitive system is governed by the system 

q = i (26a) 

z' = h{z,v) (26b) 

q = Cj(z,v,t), (26c) 

whereas the equations defining a meminductive system are 

u = p(u,i) (27a) 

if' = v (27b) 

(p = 6(u,i,t). (27c) 

Again, with Cj(z, v) = C m (z)v, 9(u,i) = L m (u)i we get the systems analyzed by Di Ventra 
et al. in [TT] . 

Mem-systems have many potential applications, not only in electronics; see e.g. (33J, [36] 
and references therein. These systems are likely to define an active field of research in the 
near future. 



4 Nodal analysis of first order circuits 

From a computational point of view, models of the form ([I]) offer some difficulties for nu- 
merical simulation. This is due to the fact that an automatic computation of the loop and 
cutset matrices B, D is difficult to perform in practice, specially in high scale integration 
circuits. For this reason, it is often preferred to introduce the node potentials e in the model, 
and describe the circuit equations using nodal analysis. This is the case in most circuit 
simulation problems, notably in SPICE and its commercial variants, which set up the circuit 
equations using Modified Nodal Analysis (MNA); cf. [121 US HS1 EH ESI SHI SZ1 HE] - 

As detailed below, the models arising from nodal analysis naturally take the form of a 
differential-algebraic equation (DAE) [U dU [231 EZ1 EH]- The main problem in the analysis 
of such differential-algebraic models is the characterization of the index, a concept which, 
roughly speaking, extends the notion of the Kronecker-Weierstrass index of a matrix pencil 
[T5] and measures the numerical difficulties faced in simulation [U [T7] . Index one and index 
two systems require specific numerical techniques, and because of this it is important to 
characterize the circuit configurations which lead to models with these indices. Therefore, 
we undertake in this Section the index analysis of nodal models of general first order circuits, 
using the tractability index framework [HJ [2U ESI HE] and extending the results discussed 
in [TJl EH ffiH HZ] to circuits with (voltage-controlled) memcapacitors, (current-controlled) 
meminductors, and hybrid memristors. Nonlinear circuits with Chua's memristors, memca- 
pacitors and meminductors define a particular case of great interest in current applications 

rig[2S[23E2JE21E2ESlES]. 

Also from an analytical point of view it is important to characterize index one circuit 
configurations. In index one systems, all the dynamic variables of the different devices con- 
tribute to the state dimension of the problem; more precisely, in index one cases the state 
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dimension (also called the order of complexity) of a first order circuit equals the sum of the 
state orders of the devices with differential order one. By contrast, in higher index problems 
which arise from so-called topologically degenerate configurations the feasible values for these 
dynamic variables are restricted by algebraic (non-differential) constraints. In a classical 
setting the topologically degenerate configurations are VC-loops (loops just defined by ca- 
pacitors and -possibly- voltage sources) and IL-cutsets (cutsets just including inductors and 
-possibly- current sources); note that V- loops and I-cutset are excluded in well-posed prob- 
lems. As a byproduct of our analysis we introduce the notion of a topologically degenerate 
configuration for circuits including memcapacitors and meminductors. 

4.1 The nodal model 

Assuming a time-invariant setting for simplicity, the nodal model of a general first order 
circuit reads as 
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The two distinct features of nodal models are the use of node potentials e and the description 
of Kirchhoff laws as Ai = 0, v = A T e, in terms of the reduced incidence matrix A. Provided 
that the circuit is connected and that a reference node has been chosen, this matrix is defined 
as A = (ciij) with 

{1 if branch j leaves node i 
— 1 if branch j enters node i 
if branch j is not incident with node i, 

for all nodes except for the reference one. In (|28|) this matrix is partitioned by columns 
according to the electrical nature of the corresponding branches. 

Several additional remarks are in order. The size of the model (|28p owes to its very 
general nature, which accommodates a great variety of devices, namely, capacitors, induc- 
tors, current-controlled and voltage-controlled resistors, g-memristors, <£>-memristors, mem- 
capacitors, meminductors, current-controlled and voltage-controlled hybrid memristors, and 
voltage and current sources. The subscripts corresponding to these devices are c, I, r, g, m, 
w, mc, ml, hm, hw, u and j, respectively. For the sake of simplicity capacitors and inductors 
are assumed to be voltage- and current-controlled, respectively (cf. the maps 7 C and 7; in 
(I28nl) and ( |28qP ), and we eliminate the branch currents i g , i w of voltage-controlled resistors 



and 92-memristors by means of the maps 7 9 and (, respectively. The map 7 r in f)28sl) defines 
the characteristic of current-controlled resistors, whereas i s (t) and v s (t) are the excitations 
in the current and voltage sources, respectively. The maps rj, u, 6, x/j, £ are those arising in 
the characteristics (|9b|) . flllcl) . fll2cp . (II 3d) and (I14cj) . except for the fact that they account 
for the whole sets of g-memristors, memcapacitors, meminductors and hybrid memristors, 
and hence need not be scalar; the same holds, of course, for yj-memristors and the map (. 

For later use, denote by R, G, C, L, M, W, C m , L m , M h and the incremental 
resistance, conductance, capacitance, inductance, memristance, memductance, memcapaci- 
tance, meminductance, hybrid memristance and hybrid memductance matrices, defined by 
the derivatives i r , i g1 i c , i h T] im , ( Vw , u Vmc , 6 iml , ^ ihm , £ flfcw , respectively. These matrices need 
not be diagonal, meaning that full coupling is allowed within each of these sets of devices. 

4.2 Topologically nondegenerate configurations and index one models 

Semiexplicit differential-algebraic equations are defined by a system of the form 

x' = f{x,y,t) (29a) 
= g{x,y,t), (29b) 

where x G W denotes the differential or dynamic variables, y G M? stands for the algebraic 
ones, / G C 1 (M r+p+1 , W), and g G C 1 (M r+p+1 , M p ). The DAE ([29]) is said to be index one 
around a given (x*,y*,t*) satisfying (129bl) if the matrix of partial derivatives g y (x* ,y* ,t*) is 
invertible. Often, this non-singularity requirement holds everywhere. This definition applies 
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in the context of the geometric, differentiation and tractability indices, and for both analyt- 
ical and numerical purposes it is important to characterize index one systems in practice; 
detailed discussions about the different index notions can be found in [U [HI [T7J [23j |37J [39] . 
The nodal model fl28|) has a semiexplicit form, the set of dynamic variables being 

Qci Qmci ^li ^Prnli Qmi Qmh Qhmi Qhwi ^wi tymci tyhmi ^Phwi (30) 

whereas the algebraic ones are 

6, ici ifnci iui ili i"inh iri imi ihmi ihw (31) 

We address below the characterization of index one configurations for (I28p . under the as- 
sumption that certain circuit matrices are positive definite; recall that a given matrix K 
is positive definite if u T Ku > for any non-vanishing real vector u, and that this notion 
expresses mathematically a strict passivity requirement on the corresponding devices. The 
proof proceeds by showing how the non- singularity of the matrix defining index one con- 
figurations can be reduced to a form already analyzed in the context of nonlinear circuits 
without memristive devices in [39J . The result stated in Theorem [T] actually motivates the 
following definition. 

Definition 11. A first order circuit is said to be topologically nondegenerate if it does not 
display either loops defined by voltage sources, capacitors and/or memcapacitors, or cutsets 
composed of current sources, inductors and/or meminductors. 

This extends a well-known notion for RLC circuits, for which topologically nondegener- 
ate configurations preclude loops defined by voltage sources and/or capacitors, and cutsets 
composed of current sources and/or inductors. Stemming from the work of Bashkow [1] in 
the classical RLC setting, this provides a way to formulate a state space model of the circuit 
dynamics by means on the notion of a proper tree. Note that in our present setting neither 
q- or <y3-memristors, nor hybrid ones, introduce topological degeneracies. 

Theorem 1. Assume that the capacitance C , the memcapacitance C m , the inductance L and 
the meminductance L m are non-singular matrices, and that the resistance R, the conductance 
G, the memristances M, and the memductances W , are positive definite. 

Then the model is index one if and only if the circuit is topologically nondegenerate. 

Proof. The matrix of partial derivatives of the right-hand side of ( |28j) w.r.t. all variables 
but time, to be denoted by F, has the form 

where the block F22 stands for the partial derivatives of the restrictions (I28mj) - f28vp with 
respect to the algebraic variables fl3TT) . The non- singularity of F22 characterizes index one 
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(32) 



configurations, and this matrix reads as 
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the block 4 U being an identity matrix whose size is defined by the number of voltage- 
controlled hybrid memristors. The non-singularity of C, C m , L, L m obviously reduces the 
problem to the characterization of the non- singularity of 
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and, by means of a Schur reduction 
that of 



the non-singularity of this matrix amounts to 



{ AvGAT + AvWAZ + ArR^AT + AmM^Al + AhmM^A^ + A^WkAl,, A c A 
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This matrix has the structure arising in [391 Theorem 5.1(1)] (cf. eq. (5.43) there); from this 
result it follows that, in the present setting, the non-singularity of this matrix relies on the 
absence of loops composed of voltage sources, capacitors and/or memcapacitors and cutsets 
defined by current sources, inductors and/or meminductors, as we aimed to show. 



□ 



4.3 Topological degeneracies: Index two 

In presence of the topologically degenerate configurations discussed above (i.e. loops defined 
by voltage sources, capacitors and/or memcapacitors, or cutsets composed of current sources, 
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inductors and/or meminductors) , the index one condition for the nodal system ( 128|) fails. 
In this situation, for both analytical and numerical purposes it is important to characterize 
whether the model is index two or not. The index two notion for a DAE is more intricate than 
the index one concept introduced above. Again, the reader is referred to [U HH [T71 [23J [37J ES] 
for different approaches to the index notion. 

In particular, the tractability index notion, together with the projector-based framework 
supported on it [HI [2U [251 ISfE ESI H21 EE], has been proved to be a valuable tool in circuit 
simulation [121 US [TU [2?1 [391 SZl SH]- In order to introduce this notion, we look at (F5B]) as 
a semilinear problem of the form 



where E is a block-diagonal matrix block-diag{/, 0}, and z joins together the variables 
denoted by x and y in (I29j) . Consider the matrix pencil XE — F [13], F being the matrix 
of partial derivatives f z . As detailed in the references above, the pencil is said to have 
tractability index one if E\ = E — FQ is a non-singular matrix. In turn, if E\ is singular, 
we let Q\ be any projector onto ker£"i, and the pencil is said to have tractability index two 
if E2 = E\ — F\Q\ is non-singular, where F\ = F(I — Q). 

Iteratively, this approach provides a general index notion which can be shown to equal 
the Kronecker-Weierstrass (or nilpotency) index of the pencil, being well-suited for compu- 
tational purposes. Moreover, this concept can be extended to nonlinear and/or time- varying 
settings under suitable assumptions on the system operators; restricting the attention to 
DAEs of the form ( |33l) in an index two context, these assumptions amount to requiring that 
Qi be a continuous projector onto the kernel of E\(z). Supported on these ideas, we show 
below that the nodal model (J28l) is indeed index two in the presence of degenerate config- 
urations, under additional passivity requirements. Note that in the index two context the 
normal tree method of Bryant [HI [6] applies in order to derive a state space equation. 

Theorem 2. Assume that the capacitance C , the memcapacitance C m , the inductance L, 
the meminductance L m , the resistance R, the conductance G, the memristances M , and 
the memductances W , are positive definite. 

Then the nodal model has tractability index two in the presence of topologically 
degenerate configurations. 

Proof. As indicated above, the model (1281) has the form depicted in (1331) with E =block- 
diag{7, 0} and F being the matrix or partial derivatives displayed in ( 1321) . Letting Q be a 
projector onto the kernel of E with the structure block-diag{0, /}, we arrive at 



Note, incidentally, that this makes it clear why the index one condition arising in Theorem [T] 
relies on the non-singularity of F22. As indicated in the proof of Theorem (U the presence of 
loops defined by voltage sources, capacitors and/or memcapacitors, and/or cutsets composed 



Ez' = f(z,t) 



(33) 
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of current sources, inductors and/or meminductors makes F 22 (and hence E\) a singular 
matrix. 

In order to construct a projector Q\ onto ker£i, we need to detail the structure of F\ 2 
and F21 in ( 132]) . These blocks are defined by the derivatives of (I28al) - fl281j) w.r.t. the algebraic 
variables (I3T]) and the derivatives of the restrictions f l28mp - fl28vp w.r.t. the dynamic variables 
f pOj) . respectively, and they read as 
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with K = A w d(/dcp w , K\ = -du/dip mc , K 2 = -d9/dq m i, K 3 
K 5 = dip/diphm, K 6 = -d£,/dq hwi K 7 = -d£/dip hw . 

Now, the projector Qi may be chosen to have the structure 



drj/dq m , K 4 = dip/dq hm , 
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where Q a and Qb are 
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being non-trivial in 



Here Q is a projector onto ker (A c A Ti 
the presence of cutsets defined by inductors, meminductors and/or current sources, whereas 

/ Qn Q12 Q13 \ 

Q — Q21 Q22 Q23 
\ Q31 Q32 Q33 J 

is a projector onto ker (A c A mc A u ), which does not vanish in the presence of loops composed 
of capacitors, memcapacitors and/or voltage sources. 

Additionally, the matrix F x = F(I — Q) has the expression 





F21 



This gives E 2 = E\ — F\Q\ the form 
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and, except for the — sign, the lower-right block (which characterizes the non-singularity of 
E 2 ) reads 
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Now, by means of a Schur reduction the non-singularity of this matrix is easily proved 
equivalent to that of 



A X M X A\ + A 2 M 2 A\Q A 
-M 3 AJ 




with A 1 = {A g A w A r A m A hm A hw ), M x = block-diag(G, W, Br 1 , M" 1 , Mr 1 , Wh), A 2 = 
(Ai A ml ), M 2 = block-diag(L- 1 ,L m 1 ), A 3 = (A c A mc ), M 3 = block-diag(C, C m ), A A = A u 



and 




Q= | Qn ? 12 | , Q 

Q21 Q22 

The form of this matrix amounts to that arising in the index-two analysis of Augmented 
Nodal Analysis models of nonlinear circuits without memristive devices, which is proved 
in [39j Theorem 5.1(2)] to be non-singular provided that Mi, M 2 , M3 are positive definite. 
In this case, the definiteness of these matrices follows from the assumption that the circuit 
matrices are positive definite and the proof is complete. 

□ 



In problems without memristive devices, Theorems [T] and [2] particularize to the results ob- 
tained in [12j [39j ST] . Our results also extend the characterization derived in [13] for cases 
including Chua's memristors. The scope of the general index characterization here presented 
covers all types of devices with differential order one, and therefore applies also to circuits 
with q- and yj-memristors, voltage-controlled memcapacitors, current-controlled meminduc- 
tors, and hybrid memristors. These results should be useful in future analytical or numerical 
studies of general first order circuits. 



5 Concluding remarks 

We have presented in this paper a comprehensive taxonomy of a variety of devices which 
have arisen in nonlinear circuit theory in the last few years. This taxonomy is organized 
around the notions of the differential and the state order of a device. An exhaustive list of 
possibly nonlinear devices with differential order one has been discussed: besides capacitors 
and inductors, these include q- and y9-memristors, memcapacitors, meminductors, and also 
the hybrid memristors here introduced. Hybrid memristors display a characteristic relating 
all four fundamental circuit variables, and account for devices in which memory effects of 
resistive, capacitive and inductive nature coexist. All these devices are discussed using fully 
nonlinear characteristics, and particularize to Chua's memristors and to the memcapacitors 
and meminductors of Di Ventra et ai, respectively, when the corresponding constitutive 
relations are linear in certain variables. A detailed analysis of the differential-algebraic 
index of circuits including all possible types of first order devices is also included. 
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Many aspects remain open and define lines for future investigation. These include nu- 
merical issues, modelling aspects involving e.g. branch-oriented systems and hybrid analysis, 
or dynamical properties related to the nature of these circuits' operating points, their stabil- 
ity, bifurcations, as well as the eventual existence and characterization of periodic solutions, 
oscillations or chaotic effects. Higher order devices and mem-systems are also in the scope 
of future research. 
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